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In the framework of nonrelativistic QCD, we study the radiative decays of bottomonia into char- 
monia, including T — * XcJl, T — » ??c7, Vb J/i'li and xu ~* J/ipf- We give predictions for their 
branching ratios with numerical calculations. E.g., we predict the branching ratio for r/t — * J/ip'y 
is about 1 x 10 -7 . As a phenomenological model study, we further extend our calculation to the 
radiative decays of bottomonia into light mesons by assuming the /2(1270), /2(1525) and other 
light mesons to be described by nonrelativistic qq (q — u, d, s) bound states with constituent quark 
masses. The calculated branching ratios for T — * /a(1270)7 and T — » / 2 (1525)7 are roughly con- 
sistent with the CLEO data. Comparisons with radiative decays of charmonium into light mesons 
such as J/ip — > /2(1270)7 are also given. In all calculations the QED contributions are taken into 
account and found to be significant in some processes. 

! PACS numbers: 12.38.Bx; 13.25.Hw; 14.40.Gx 
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CN '. I. INTRODUCTION 

Radiative decays of bottomonmm (e.g. T, 77b, Xbj) into charmonium are expected to be described by nonrelativistic 
I quantum chromodynamics (NRQCD), since both bottomonium and charmonium are made of heavy quark and heavy 
antiquark, and are nonrelativistic bound states. For heavy quarkonium decay and production, the rates can be 
factorized into a short-distance part, which can be calculated in QCD perturbatively, and a long-distance part, which 
are governed by nonperturbative QCD dynamics [l[. Therefore, radiative decays of bottomonium into charmonium 
may provide a useful test for NRQCD factorization, which is assumed to hold also for these specific exclusive processes, 
and may also provide some practical estimates for decays such as rjb — ► J/ifij, which might be useful in search for the 
not yet discovered rjb meson. As a phenomenological model study, we further extend our calculation to the radiative 
decays of bottomonia into light mesons by assuming the /2(1270), f 2^525) and other light mesons to be described by 
^p-i. nonrelativistic qq (q = u, d, s) bound states with constituent quark masses m q '(q — u, d) — 350 MeV, m s = 500 MeV 
O ,' as in constituent quark models. These radiative decays are known as the gluon rich channels, and regarded as a good 
place to investigate the interactions between quarks and gluons in these OZI forbidden processes, and there have been 
some earlier work discussing these processes (see, e.g.,0, Q). In this paper, as our previous work we will perform a 
, complete numerical calculation for the quark-gluon loop diagrams involved in these processes, and we will also include 
contributions from QED diagrams in the same processes. 

We adopt the assumption that both heavy quarkonium and light mesons are described by the color-singlet non- 
relativistic wave functions. Based on this assumption, we study T — > XcJ7, Y — ► ?7c7, T — > fjj, T — > 777, J/ip — > /,/7, 
J/ip -> ?77, xbj -> J/ip(p,uJ,(f>)-f and -q b -> J/ip(p,ui,4>)j etc. 

The rest of this paper is as follows. In section II, we will give the descriptions and main techniques in our calculations, 
and then make predictions for the decay rates of T — > XcJl-, T - > ?7c7, Vb ~^ J/^l, and Xbj — * J 1^1- Then, in the 
following section, we will generalize this method to those processes in which the final states are light mesons. Finally, 
we will summary all the results in section IV. 
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II. BOTTOMONIUM RADIATIVE DECAYS TO CHARMONIUM 



In this section, we will study the radiative decays of bottomonium into charmonium. In NRQCD, heavy quarkonium 
wave function is described by a Fock state expansion in terms of the relative velocity v between the quark and 
antiquark, and the leading term is a color-singlet QQ state , which has the same quantum numbers as the physical 
heavy quarkonium. In certain processes, the non-leading terms with color-octet QQ pair and soft gluons may make 
dominant contributions. E.g., in the T radiative decays to light quark jets T — > qqj the color-octet contribution could 
be larger that the color-singlet contribution (depending on the estimates of the color-octet matrix elements) [5[ . In the 
radiative decays of bottomonium into charmonium, the short distance transitions of a color-octet bb into a color-octet 
cc by emitting a photon are shown in Fig.l, where q = c, and qq are in color-octet ( 3 <Si) 8 or ( 3 Pj) 8 . Compared with 
the case of T radiative decays to light quark jets T — ► qqj (see Ref.[j| for an estimate of the color-octet contributions), 




(b) (c) 
FIG. 1: Feynman diagrams for transitions from color octet bb to color-octet qq (where q — c) by emitting a photon 





FIG. 2: Typical QCD Feynman diagrams for heavy quarkonium ( 3 Si) radiative decays into { 3 Pj) mesons 

here the contribution of color-octet cc is greatly suppressed by the smallness of the color-octet matrix elements of 
( 3 Si) 8 or ( 3 Pj) 8 (note that the color-octet matrix element of ( 3 Si) 8 is only 1% of that of color-singlet ( 3 S'i) 1 for J/ip. 
Therefore, we will neglect the color-octet contributions, and only concentrate on the color-singlet description of heavy 
quarkonia in the following calculations. 

A. General results 

In the nonrelativistic approximation, the T radiative decay into a color-singlet cc pair, which subsequently 
hadronizes into charmonium, can be described by the diagrams in Fig. 2, and the amplitude can be expressed as[|| 
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where (3Z;3« | 1)= 5u/^fW~ c , (si;S2 | S'S'z) , and (LL Z ;SS Z \ JJ Z ) are respectively the color-SU(3), spin-SU(2), and 
angular momentum Clebsch-Gordan coefficients for QQ pairs projecting on appropriate bound states. A(bj(pb) + 
bk(pb) — > Qi(Pc) + Qi(Pc)) is the decay amplitude for QQ production and ^4° is the derivative of the amplitude with 
respect to the relative momentum between the quark and anti-quark in the bound state. The coefficients Ql t and Cl 
can be related to the radial wave function of the bound states or its derivative with respect to the relative spacing as 



C S = ±\R S (0)\ 2 , C P = ^\R' p (0)\ 2 . 



(2) 



The spin projection operators Pss z (Pi l) which describe production of quarkonium are expressed in terms of quark 



FIG. 3: Typical QED Feynman diagrams for heavy quarkonium ( 3 Si) radiative decays into ( 3 Pj) mesons 



and anti-quark spinors as@, Q: 

PssAP^)=^2 v (^-q^2)u(^+q,s 1 )(s 1 ;s 2 \SS z ), (3) 

We list the spin projection operators and their derivatives with respect to the relative momentum, which will be used 
in the calculations, as 

P 0O (p,0) = JL 75 (^+2m), (4) 
PlSzM = f(S,)V + 2m), (5) 

Pf s >,0) = -^-[l a f(S,)U + 2m)-y-2m) f(S z ) 7 a }. (6) 

And the spin projection operators which describe the annihilation of quarkonium are the complex conjugate of the 
corresponding operators for production. The polarization vectors for the 3 Pj states are shown below: 

^e*"(L z )e*P(S z ){lL z -,lS z \lJ z ) = ~^^7* m , (7) 

LzSz 

J2e*%*y (S z ){lL z -lS z \OO) = -L^ + ^ft, (8) 

LzSz * 6 

Y / £* a (Lz)e* p (S z )(lL z ;lS z \2J z ) = e* a ?{J z ), (9) 

LzSz 

where p is the momentum of P-wave quarkonium, and M is the mass of the corresponding quarkonium. e\(J z ) are 
the polarization vectors for J = 1. e a @(J z ) are the polarization vectors for J = 2, which are symmetric under the 
exchange a «-►/?. 

The QCD Feynman diagrams of T — > 7?7c(Xcj)> i n which the cc are produced through gluons, are shown in FIG[2l 
while the QED Feynman diagrams, in which the cc are produced through the photon, are shown in FIGE1 In the 
calculation, we use FeynCalc [f| for the tensor reduction and LoopTools[9j for the numerical evaluation of infrared 
safe integrals. We follow the way in Ref.[10( to deal with five-point functions and high tensor loop integrals that can 
not be calculated by LoopTools and FeynCalc, such as 

p f jDi kgkpkpkg 

afip ° ~ J fc 2 [(* + p c ) 2 - ml] (k + 2p c n(p b + kf - ml] [(k + 2p c -p b y-mlY 1 Uj 

where p c is the momentum of c quark, p b the momentum of b quark, m c the charm quark mass, and m b the bottom 
quark mass. 

In the numerical calculations, the quark masses are taken to be m b = 4/7 GeV, m c = 1.5 GeV, the wave functions at 
the origin can be found from potential model calculations in Ref.p]]: |7^ b (0)| 2 = 6.477 GeV 3 , \Tlf (0)| 2 = 0.81 GeV 3 , 
l^p b (°)| 2 = !- 417 GeV 5 , |ft' cE P (0)| 2 = 0.075 GeV 5 . In the bottomonium decay, the strong coupling constant is chosen 
as a s (2m b ) — 0.19. The numerical results of radiative decays of bottomonium into charmonium are listed in Table I. 
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process 


T -> Xc27 




T -> Xc07 


T -> 77c7 


BRqcd 
BRqcd+qed 


5.1 x 10" 6 
5.6 x 10" 6 


4.5 x 10" 6 
9.8 x 10" 6 


4.0 x 10" 6 
3.2 x 10" 6 


2.9 x 10" 5 
4.9 x 10" 5 


process 


Xb2 -> 


Xbi -> J/V>7 


Xw —> J/ipl 


1% -> J/V>7 


r QCD (GeV) 
Fqcd+qed (GeV) 
Tq-B-D (GeV) 


2.7 x 10" 10 
3.6 x 10 -10 

3.8 x 10" 11 


3.8 x 10" 10 
3.7 x 10~ 10 
3.3 x 10" 12 


5.0 x 10" 10 
1.3 x 10 -10 
1.3 x 10" 10 


2.8 x 10" 9 
9.6 x 10 -10 
1.2 x 10" 9 



TABLE I: Decay widths and branching ratios for radiative decays of bottomonium into charmonium. The decay widths V are 
in units of GeV, and the branching ratios BR are given for the T. 



rib radiative decay to J/tp 



The rib meson is the only one among the low lying bottomonium states that has not been observed experimentally. 
To search for the rjb meson a number of decay channels have been suggested, e.g., decays into the J/ipJ/ip and 

DD [Sill fl7j |. In any case, the radiative decay rjb — ► J/ipj should be a useful channel for the 17& in view of the 
cleanness of the signal (this possibility has also been considered in Ref . fl8|| ) . 

From Table I, we can see that for the i]b decay 1% — > J/tp'y the QCD and QED contributions are comparable but 
destructive, and, as a result, the decay width of 1% — > J/ipj is only 9.6 x 10~ 10 GeV. In order to know the branching 
ratio of this decay channel, we should have an estimate for the rjb total width. In fact, we can estimate its total width 
through Ttotivb) ~ T(?7& — * gg) [J- For T(rib — > gg), with next to leading order (NLO) QCD radiative corrections, we 
have 



r(% -> gg) 



\n s (0)\ 2 C F a 2 s (2mb) 
2m\ 



1 
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~24~ 
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(11) 



With the parameters used above, we can get TtotiVb) ~ H-4 MeV. Then the branching ratio is Br(i]b — > jJ/ip) ~ 
8.4 x 10~ 8 . If we use the leading order formula in Eq. p4|) . the decay width is TtotiVb) ~ 7.1 MeV, and the branching 
ratio becomes Br(r)b — ► jJ/ip) ~ 1.4 x 10~ 7 . 

On the other hand, with the spin symmetry in the nonrelativistic limit (v = 0), the rjb wave function at the origin 
|7?. s (0)| 2 can be determined from the T leptonic width, 



r(T^e+e~) = N c Q 2 b a 



l^(Q)l 2 

3to? 



and the rjb total width is then related to the T leptonic width, 



Ttot(Vb) 



3C F a 2 (2m b 



1 



! 199 _ 13tt 2 
I 18 24 



1 - 



C A 



"3T 



2N c Q 2 b a 2 



1 



16a s 
3tt 



-r(T 



■)• 



(12) 



(13) 



Using m b = 4.7 GeV, a s (2m b ) = 0.19, and experimental data T(T -> e+e~~) = 1.340 ± 0.018 KeVjlJ], we can get 
FtotiVb) ~ 13.0 MeV. Then the branching ratio is Br(r]b — ► jJ/ip) ~ 7.3 x 10~ 8 . If we use the leading order formula 
in Ea.([i"i"1) and Eq.©, the T to t{r)b) « 5.45 MeV, the branching ratio is Br(ri b -> 7J/VO « 1.7 x 10" 7 . 

In any case, we find that the branching ratio Br(rib — > jJ/tp) is of order 1 x 10 -7 . This small number makes it 
quite difficult to search for r/ b through this decay channel. 



C. Helicity ratios with Xci and Xc2 

We give predictions for branching ratios for different helicity states in T — > XcJl decays. As in Ref. we choose 
the moving direction of \cj as the z-axis, and introduce three polarization vectors: 

= -±(0,M,0), 
c^(-l) = -^(0,l,-i,0) 

w»(0) = — (|k|,0,0, (14) 
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QCD 


QCD+QED 




QCD 


QCD+QED 


y 2 (T -> 7x^2) 


0.37 
0.14 


0.38 
0.14 


^ 2 (T -» 7Xd) 


0.064 


0.075 



TABLE II: Results for T -+ XcJj(J = 1, 2) with different helicity states 



thus we can characterize the tensor e MI/ (A) of %c2 

e Q/3 (2)=w Q (lV /3 (l) 

e^(l) = -^(w Q (l)^(0)+cj a (0)a; /3 (l)) 
v2 

e Q/3 (0) = ^K(-l)^(l) + 2^(0)^(0) +w tt (l)^(-l)) 
V6 

e Q/3 (~l) = ^(w Q (0)^(-l)+w Q (-l)w' 3 (0)) 
v2 

e af) {-2)=u a {-l)uP{-l) (15) 
The helicity ratios are introduced as 

, 2 = M and ^ _ M 2 



aol 2 '' M 2 ' 



where ax, A = 0, 1, 2, are the normalized helicity amplitudes, which satisfy |ao| 2 + |ai| 2 + |a,2 1 2 = 1. Namely, |<z,\| 2 
is the probability of the final state meson with helicity ±A. Then the ratios x 2 and y 2 only depend on the mass ratio 
rric/mb- With the same choice of parameters, we predict ratios for different helicities in Table II. 



III. RADIATIVE DECAYS OF HEAVY QUARKONIUM INTO LIGHT MESONS 



As a purely phenomenological model-dependent study, In this section we will extend our calculations performed 
above for radiative decays of bottomonia into charmonia to the radiative decays of bottomonia into light mesons. Our 
assumption is that the light mesons such as the /2(1270), /2(1525), and /i(1285) can be described by nonrelativistic 
qq (q = u,d,s) bound states with constituent quark masses. 

In the numerical calculations, the light quark masses are taken to be m s = 0.50 GeV, m u — — 0.35 GeV. The 
parameters for the heavy quarks are the same as that used in section II,|7£g (0)| 2 = 6.477 GeV 3 , |7?.^ c (0)| 2 = 0.81 GeV 3 , 
|ft$ 5 (0)| 2 = 1.417 GeV 5 , |7^ E (0)| 2 = 0.075 GeV 5 , m b = 4.7 GeV, and m c = 1.5 GeV. The strong coupling constant 
is chosen as a s = 0.19 and a s — 0.26 in bottomonium and charmonium decays respectively. 

As widely accepted assignments we assume that / 2 (1270) and /i(1285) are mainly composed of (uu + dd)/y/2 
(neglecting the mixing with ss for simplicity). But for /o(980), there are many possible assignments such as the 
tetraquark state, the KK molecule, and the P-wave ss dominated state (for related discussions on /o(980) and other 
scalar mesons, see, e.g., the topical review-note on scalar mesons in [l4j and Since experimental data show that 

Df — > /o(980)7r + has a large branching ratio (BR)[l3|, here we assign /o(980) as an ss dominated P-wave state as a 
tentative choice (we do not try to justify this assignment). 

As to the wave functions at the origin of light mesons, it is very difficult to determine them without any doubt. 
Using the theoretical expression for the widths of fi — ► 77 [l|, 



r (th) 6AT ^JM^Ln 8as ^ 

r / 2 (1270)-, 77 = —(Qu + Q d ) " m 4 i 1 ~ 3^) 

r (th) 12V C 4 2 \K> P (0)\ 2 H 8a s 2 

^(1525)^ = —Qs& m 4 ' ( 16 ) 

where N c =3 is the color number, a = 1/137, and fitting them with their experimental values 2.6 KeV and 0.081 KeV 
for / 2 (1270) and /a (1525) respectively [3 , we get 

1-3 /-t„Tr5 



\KT{0)\ = 1-58 x 10~ 3 GeV 

) /s. 

-p 



|^ 5 (0)| 2 = 2.23 x 10" 3 GeV 5 . (17) 
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process 


T - * 7/0 (tot) 


T — > 7/1 (1285) 


T -> 7/2(1270) 


T -> 7/0 (980) 


7^ (1420) 


7/2(1525) 


BR?r 


7.2 x 10" 5 


2.6 x 10" 5 


7.1 x 10" 5 


2.2 x 10" 5 


1.0 x 10" 5 


2.2 x 10" 5 


ggQCD+QED 


6.4 x 10" 5 


3.9 x 10" 5 


6.3 x 10 -5 


2.0 x 10" 5 


1.2 x 10" 5 


2.0 x 10" 5 


BR ex 


\ 


\ 


1.00 ±0.10 x 10" 4 


< 3 x 10" 5 


\ 


3.7t\i x 10 -5 


process 


J ftp -> 7/0 (nii) 


J/i>^ 7/1 (1285) 


J/^ 7/2 (1270) 


J/i> -> 7/o (980) 


J/V^ 7/1(1420) 


J/i> -> 7/2(1525) 


BR? h CD 


2.0 x 1(T 3 


2.0 x 10" 3 


2.5 x 10" 3 


6.7 x 10" 4 


7.5 x 10~ 4 


8.5 x 10~ 4 


BR QCD+QED 


1.8 x 10~ 3 


2.7 x 10~ 3 


2.4 x 10 -3 


6.5 x 10~ 4 


8.5 x 10~ 4 


8.5 x 10~ 4 


BR ex 


\ 


(6.1 ±0.8) x 10~ 4 


(13.8 ± 1.4) x 10" 4 


\ 


(7.9 ±1.3) x 10" 4 


(4.5l°;I) x 10- 4 



TABLE III: Numerical results for T(J/ip) -> 7/j . 



If we use the leading order formula in Eq (|16p . then 

|^"(0)| 2 = 6.6 x 1CT 4 GeV 5 

\K'f(0)\ 2 = 1.1 x 1(T 3 GeV 5 . (18) 

For the vector mesons, the wave functions at the origin may be determined from their leptonic decay V — » 
e+e - (V — <j),p) widths. Using 

r(#1020) -► e+e") = V C Q^ 2 ^P(1 _ ^£)2 = (1 27 ± .04) KeV, (19) 

3m, 37r 



we can get 



|7^ s (0)| 2 = 0.11 GeV 3 , 

\K S S S (0)\ 2 = 0.19 GeV 3 . (20) 



If we use the leading order formula in Eq . (fl9|) . then 



|7^"(G)| 2 = 0.032 GeV 3 , 

\TI S S S (0)\ 2 = 0.054 GeV 3 . (21) 

The wave functions at the origin of light mesons can also be determined from potential models [20] ] . From experi- 
mental data, AE = M(2S) - M(1S) is 675 MeV, 638 MeV, 661 MeV, 589 MeV, and 563 MeV for p, u>, (f>, J/tp, and 
T respectively. In the logarithmic potential, AE is independent of quark masses. So we may select the logarithmic 
potential, which gives 

\R S (0)\ 2 cx m 3 / 2 

\R P (0)\ 2 oc m 5 / 2 (22) 
With \nf (0)| 2 = 0.81 GeV 3 and \ll'f{0)\ 2 = 0.075 GeV 5 1 . Then we can get 



3/2 

|7^ fi (0)| 2 = ( — ) |ftf(0)| 2 = 0.091GeV 3 



m, 



3/2 

\TZ S S S (0)\ 2 = (— ) |ftf(0)| 2 =0.156GeV 3 



ni, 



.5/2 

|7^ fi (0)| 2 = ( — ) |7^(0)| 2 = 1.97 x 10- 3 GeV 5 
1 m r ' 



f \ 5 ^ 2 

Tl'Z s (0)\ 2 = — ) |^ e (0)| 2 = 4.81 x 10- 3 GeV 5 (23) 
\m c J 



1 The wavefunction at the origin in the logarithmic potential for cc is |7£= c (0)| 2 = 0.815 GeV 3 and |7^ c (0)| 2 = 0.078 GeV 5 . It is 
consistent with the B-T potential result that was used here. 
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process 


Xbi -> 7P 


Xbi -» 7P 


XbO -> 7P 


7?6 — > 7P 


r Q ci3(GeV) 
Tqcd+qed (GeV) 


1.1 x 10" 10 
6.8 x 10" 10 


2.2 x 10" 10 
2.2 x 10" 10 


3.5 x lO -11 
8.4 x 10" 10 


1.2 x 10" 10 
1.1 x 10" 8 


process 


Xbi —> 7^ 


Xm -> 7^ 


Xbo -> 7^ 


?7b — ► 7W 


r Q c.D(GeV) 
Fqcd+qed (GeV) 


1.2 x 10" 11 
7.6 x 10" 11 


2.4 x 10" 11 
2.4 x 10" 11 


3.8 x 10" 12 
9.3 x 10" 11 


1.3 x 10" 11 
1.2 x 10" 9 


process 


Xb2 -» 70 


Xbi -> 70 


Xbo -» 70 


Vb -» 70 


r Q c.D(GeV) 
Fqcd+qed (GeV) 


3.6 x 10 -11 
1.3 x 10" 10 


5.8 x 10" 11 
5.8 x 10" 11 


1.6 x 10" 11 
7.5 x 10 -11 


5.9 x 10" 11 
1.2 x 10" 9 


process 


Xc2 -> 7P 


Xd -» 7P 


Xco -> 7P 


J? c -> 7P 


BR? h CD 

gpQCD+QED 


1.3 x 10" 5 
3.8 x 10" 5 


4.1 x 10" 5 

4.2 x 10" 5 


3.2 x 10" 6 
2.0 x 10" 6 


1.5 x 10" 6 
4.5 x 10" 6 


process 




Xci -» 7^ 


Xco -» 7^ 


r/ c — > 70; 


BRfr 

BR QCD+QED 


1.5 x 10" 6 
4.2 x 1(T 6 


4.6 x 10" 6 

4.7 x 10 -6 


3.5 x 10" 7 
2.2 x 10 -7 


1.7 x 10" 7 
5.0 x 10" 7 


process 


Xc2 —> 70 


Xci -» 70 


Xco -» 70 


Vc -> 70 


BR? h CD 


3.3 x 1(T 6 
6.5 x 10" 6 


1.1 x 10" 5 
1.1 x 10" 5 


1.3 x 10" 6 
3.0 x 10" 8 


7.1 x 10" 7 
4.0 x 10" 7 



TABLE IV: Predicted decay widths for Xbjiflb) — > 7p(w, 0). 



m„(GeV) 




QCD 


QCD+QED 




QCD 


QCD+QED 




QCD 


QCD+QED 


0.35 


* 2 (T -> 7/2) 


0.023 


0.023 


y 2 (T^ 7 / 2 ) 


0.0013 


0.0014 


* 2 (T - 7/O 


0.00069 


0.00010 


0.635 


z 2 (T - 7/2) 


0.073 


0.075 


y 2 (T~, 7 / 2 ) 


0.0095 


0.010 


x 2 (T -» 7/1) 


0.0047 


0.0067 


0.50 


* 2 (T - 7/2) 


0.045 


0.046 


y 2 (T- 7 / 2 ) 


0.0043 


0.0045 


* 2 (T - 7/1) 


0.0018 


0.0023 


0.763 


^ 2 (T - 7/2) 


0.10 


0.11 


y 2 (T- 7 / 2 ) 


0.017 


0.018 


* 2 (T - 7/0 


0.0087 


0.010 


0.35 


^VA/> - 7/2) 


0.21 


0.21 


y 2 (J/^ - 7/2) 


0.055 


0.057 


x 2 ( J/i> - 7/x) 


0.026 


0.029 


0.635 


z 2 (J/V> - 7/2) 


0.62 


0.62 


y 2 (J/^^7/2) 


0.33 


0.33 


x 2 (J/^ - 7/O 


0.13 


0.14 


0.50 




0.40 


0.40 


y 2 (J/^- >7 / 2 ) 


0.16 


0.16 


x 2 (J/tA - 7/O 


0.072 


0.074 


0.763 


x 2 (J/^ -> 7 ys) 


0.86 


0.86 


yV/^7/2) 


0.57 


0.57 


z 2 (J/V^ 7 /0 


0.21 


0.22 



TABLE V: Results for T — ► XeJ 7 with different helicity states 



In the numerical calculation, the parameters are taken to be \n s s s {0)\ 2 = 0.054 GeV 3 , \TZ r ^{0)\ 2 = 0.032 GeV 3 , 
|7^(0)| 2 = 1.1 x 10" 3 GeV 5 , |7?.£™(0)| 2 = 6.6 x 10~ 4 GeV 5 . The numerical results are shown in Table III and Table 
IV. 

The branching ratio of T radiative decay into a light meson is smaller than the corresponding branching ratio of 
J/t/j by a factor of 

BR(T -> 7 M) (Q b \ 2 (m c \ 2 a(2m b ) 
BR(J/tjj —¥ 7 M) ~ \Qc) \mj a(2m c ) ~ ' 1 ' 

We can find this theoretical ratio is 0.013 ~ 0.036. The experimental ratio is 0.072 for / 2 (1270), and 0.082 for 
/2(1525). The corresponding theoretical ratios are 0.026 and 0.024 for /a(1270) and /2(1525) respectively. If we use 
larger constituent quark masses, e.g. m u = m<j = M(1270)/2, m s — M(1525)/2, the ratios are 0.018 and 0.018 for 
/ 2 (1270) and $(1525) respectively. 

With the same parameters, we give the branching ratios for different helicity states in Table V. The corresponding 
values of helicity parameters x and y are x = 0.46, y = 0.23. Recently, new experimental data for the contributions of 
different helicities in process J/tp — > 7/ 2(1270) have been given by the BES Collaboration [? ]: x = 0.89 ± 0.02 ± 0.10 
and y = 0.46 ± 0.02 ± 0.17 (see alsoflj]). It is about 2 times larger than our results. But if we use a larger constituent 
quark mass, e.g. m u — M(1270)/2, we will get substantially increased values x = 0.79 and y = 0.58 (also see Ref.0). 
We emphasize that the helicity parameters are very sensitive to the light quark masses, and hence very useful in 
clarifying the decay mechanisms. Note that if m u /m c — > 0, we will have x — > and y — > 0, but this is inconsistent 
with data. 
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IV. SUMMARY 



In this paper, we mainly investigate the radiative decays of bottomonium into charmonium, such as T — > XcJli 
Xbj — > J/ipJ, Vb —> J 1^1 an d T — > ry c 7 based on the NRQCD approach. Based on our numerical calculations, 
we predict that the branching ratios of T — > XcJl an d decay widths of Xbj ~ * J I^Pl- All the above processes are 
perturbative calculable, and it is a good way to test NRQCD. 

We next focus on the cases of heavy quarkonium radiative decays into light mesons, including T(J/ip) — > fjj and 
Xbj -> p(u, (f>)-f et ac. 

In this work, we also find that the QED effects in some radiative processes are really significant. For T — > 7Xcj 
decay, the pure electromagnetic process only affects the final results for J = 0, 2 a little, but for the J = 1 state the 
result may change by a factor of 2. The same results will be seen in the decays of T — > jfj and J/ip — ► 7/,/. As the 
cases of XbJ decays, especially for the process Xbj </>)7, QED process may give dominant contributions. 
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